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\f} \ With the consideration of spherical symmetry for the potential and 

mass function, one-dimensional solutions of non-relativistic Schròdinger 
equations with spatially varying effective mass are successfully extended 

T _— i [ to arbitrary dimensions within the frame of recently developed elegant non- 

perturbative technique, where the BenDaniel-Duke effective Hamiltonian 
in one-dimension is assumed like the unperturbed piece, leading to well- 

""qJ known solutions, whereas the modification term due to possible use of other 

i ' effective Hamiltonians in one-dimension and, together with, the corrections 

coming from the treatments in higher dimensions are considered as an 
additional term like the perturbation. Application of the model and its 
qh generalization for the completeness are discussed. 
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1 Introduction 



Gaining confidence from the successful applications [1] of the recently developed 
simple model [2] in different fieids of physics, we have investigated [3] the relation 
between the solutions of physically acceptable effective mass Hamiltonians pro- 
posed in the literature for the treatment of one dimensional problems. Using the 
spirit of the prescription suggested in [3] , we aim here to tackle the more difficuit 
problem of generating exact solutions for position-dependent mass Schròdinger 
equations (PDMSE) in N— dimension, as the most of the related works in the 
literature have been devoted to one-dimensional systems except the ones in [4]. 

The concept of PDMSE is known to play an important role in different brandi 
of physics. This formalism has been extensively used in nuclei, quantum liquids, 
3 He and metal clusters. Another area wherein the such concepts provide very 
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useful tool is the study of electrònic properties of many condensed-matter sys- 
tems, such as semiconductors and quantum dots. In particular, recent progress in 
crystal-growth techniques for producing non-uniform semiconductor specimens, 
wherein the carrier effective mass depends on position, has considerably enhanced 
the interest in the theoretical description of semiconductor heterostructures. It 
has also recently been signalled in the rapidly growing field of PT-symmetric or 
more generally pseudo-Hermitian quantum mechanics. For an excel·lent recent 
review, leading to the related references, the reader is referred to [4]. 

In Section 2, the systematic treatment of N— dimensional PDMSE is pre- 
sented and closed expressions corresponding to the full wave function and energy 
spectrum for exactly solvable potentials are given. Section 3 contains the ap- 
plication of the model while the generalization of the formalism is discussed in 
Section 4. Concluding remarks are given in the last section. 



2 Theoretical Consideration 

Tracking down solvable potentials in PDMSE has always aroused interest. Apart 
from being useful in understanding of many physical phenomena, the importance 
of searching for them also stems from the fact that they very often provide a 
good starting point for undertaking perturbative calculations of more complex 
systems. 

As is well known (see, e.g., [3, 5]), the general form of radial PDMSE with 
Hermitian Hamiltonians in one-dimension gives rise to 



+ V eff {z)<5>(z) = A$(s) , (1) 



_d_ [ 1 d$(z) 
dz [M(z) dz 
where the effective potential 

(a + >y)M" M' 2 
V"{z) = V (z) + U ai (z) = Vo(z) - + («7 + ol + 7 )^r , (2) 

depends on the mass term and ambiguity parameters. Here a prime denotes 
derivative with respect to the variable, M(z) is the dimensionless form of the mass 
function m(z) = m M(z) and we have set H = 2m = 1. The effective potential is 
the sum of the real potential profile V (z) and the modification U aj (z) emerged 
from the location dependence of the effective mass. A different Hamiltonian 
leads to a different modification term. Some of them are the ones proposed by [6] 
BenDaniel-Duke (a — 7 = 0), Bastard (a = —1), Zhu-Kroemer (et = 7 = — 1/2) 
and Li-Kuhn (7 = -1/2, a = 0). 

Considering the works in [3, 7], the radial piece of PDMSE in arbitrary di- 
mensions for spherically symmetric potentials and mass functions reads 

/ rf 2 M'ÍN-l d\ L(L + N + 2) + (N-l)(N-3)/A 

[dr 2 ~M \ 2r drj V 2 + M[E-V eff 

(3) 
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where we assume that \I/(r) = F(r)G(r) which leads to 



1 ÍF" G" F'G'\ M' ÍF' G'\ 
^ l T + ~G +2 FG )-W (f + g) =U <"- K 



M 



(4) 



The effective potential in higher dimensions (N > 1) now is transformed to the 
form 

/x M'(N-l) L(L + N-2) + (N-l)(N-3)/A 



M 2 2r 



Mr 2 



(5) 

in which L is the angular momentum. As the one-dimensional calculations require 
N — 1 and L = 0, Eq.(5) reduces in this case to U e ff(r) = V e ff(r) = V + U aj as 
in Ref. [3], which provides us a reliable testing ground. 

Keeping in mind the spirit of the technique used simply in [3], we split Eq. 
(4) in two parts deviating from the treatments in [4] 



W 2 {r) - 



W{r) 



V (r)-e, W 



F' 



MF 



(6) 



where e is the corresponding energy of the required quantum state F n (n = 
0,1,2,...) for Vo which is assumed in this model as an exactly solvable mass- 
dependent potential, and 



AW 2 (r)- 
where 



AW(r) 



M 



AV(r) = Ua~f(r) 



+2W(r)AW(r) = AV(r)-AE, AW(r) 



M' (N 



G> 



MG 



, (7) 



1) L(L + iV-2) + (iV-l)(Ar-3)/4 
M 2 2r Mr 2 [ ' 

Note that the total energy appearing in (4) is E = e + AE and, in one-dimension 
the modification term AV becomes U aj as in Ref. [3]. This clarifies that the 
corrections due to the higher dimensions arise because of the second and third 
term on RHS of Eq. (8). 

From the present theoretical consideration, Eq. (6) has an algebraic solution 
leading to closed analytical expressions for the wave functions and energy eigen- 
values, hence one needs to solve Eq. (7) exactly. To proceed further, with the 
consideration of relativistic Dirac equations having no ambiguity parameters, we 
confidently choose 

(a + 7) M' {N + 2L-Ï) 



AW(r) 



M 3 / 2 



2VMr 



(9) 



in which the second term disappears for iV = 1 as in [3] . Within the frame of Eq. 
(7), this choice leads us 



W{r)AW{r) 



M' 
2rM 2 



(a + 7)(JV- 1) 



+ (a + 7 + l)L 



AE 



(10) 
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that is the main result of the present Letter. 

From the definition of the effective potential in Eq. (2), we also note that the 
use of Eqs. (7) and (8) naturally restricts the choice of some ambiguity parameters 
yielding different physically acceptable effective mass Hamiltonians, allowing only 
a = 7 = (Ben-Daniel Duke Hamiltonian) and a = 7 = — 1/2 (Zhu-Kroemer 
Hamiltonian) cases. This observation clarifies that the unperturbed part (V ) of 
the effective potential in (6) should corresponds to the case a = 7 = 0, having 
well known solutions in one dimension, while a = 7 = — 1/2 is used to calculate 
U ai in (8). Obviously, all the corrections coming from the higher dimensions to 
the energy and well-behaved wave function terms can be systematically calculated 
for a given M with the consideration of Eqs. (7-10) in the light of corresponding 
W in (6). 

3 Application 

Recently, some researches have been devoted to the analysis of the classification of 
quantum systems with position-dependent mass regarding their exact solvability 
[3-5, and the references therein]. On a similar basis, Plastino and his co-workers 
[8] applied an approach within the supersymmetric quantum mechanical frame- 
work, for the case a = 7 = 0, to such systems and succeed to show that some 
one-dimensional systems with non-constant mass have a supersymmetric partner 
with the same effective mass. They were also able to solve exactly some particu- 
lar cases by constructing the superpotential [W(r)] from the form of the effective 
mass [M(r)] and generalize the concept of the shape invariance for these systems. 

For illustration, the superpotential expressions given by [8] for the systems 
having harmònic oscillator and Morse-like spectra can be easily used in Eq. (6) 
to serve explicit expressions for the corrections to the one-dimensional solutions 
obtained by considering the Ben-Daniel-Duke effective Hamiltonian in their [8] 
calculations. This simple investigation enables us testing our results, because all 
the corrections should disappear in case N — 1 and a = 7 = leading to the 
expressions in [8]. For clarity, this section involves only the application on the 
harmònic oscillator system. However, the generalization of the present model 
yielding self-consistent calculations, reproducing W(r) term within the model for 
any system of interest, will be discussed in the next section. 

According to Ref. [8], W(r) term in Eqs. (6) and (10) is 



W(r) = U - f sfM(z)d, 




for the system having harmònic oscillator spectra. Hence, use of Eqs. (7) through 
(10) gives 
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AE 



M' 

7mï 



{a + -y){N-l) 



+ L(a + 7 + 1) 



(N + 2L- l)u 



(-^y («+7) 



üj ( ^M(z)dz 




M 

which is the explícit form of the energy corrections for a given smooth mass. 
Clearly, it can be seen that for a constant mass M — > 1, Eq. 12 reduces to 
(n + 2L — 1)oü/2 for arbitrary dimensions [7] while in one dimension it goes to 
zero for a non-constant mass in case a = 7 = [8]. Furthermore, from Eqs. (7) 
and (9), the modification term for the corresponding wave function is 



G{r) = exp 



M(z)AW(z)dz) = r (^+2L-l)/2 M -(a +7 )/2_ 



(13) 



As Eq. (6) is analytically solvable having a closed expression for W(r) given by 
Eq. (11) reproducing explicit expressions for e and F, the corresponding total 
energy and wave function can easily be calculated through E = e + AE and 
^ = FG for the system of interest with a location dependent mass. At this 
stage it is also note that the formalism suggested here seems superior to the 
usual treatment in supersymmetric quantum theory that in principle start with 
the ground state and builds up excited state wave functions by the use of some 
linear operators (A^) whereas there is no such restriction in the present theory 
providing flexible investigations. 



4 Discussion 

Although the procedure used in the formalism seems reasonable, the use of other 
works as in the previous section for an appropriate W(r) term to solve Eq. (6) 
may be seen as a drawback of the model. To remove this seeming deficiency, we 
propose here a unified treatment within the model considering the recent work 
in [9]. 

Many of the special functions H(g) of mathematics represent solutions to 
differential equations of the form 

*JM + Q{9 f2M + R{9)H{9) ^ , (14) 

where the functions Q(g) and R(g) are well defined for any particular function 
[10]. Since in this Letter we are interested in bound state wave functions, we 
should restrict ourselves to polynomial solutions of Eq. (14). Bearing in mind 
Eq. (14), the substitution of $(z) = H[g(z)]f(z) in Eq. (1) leads to the second- 
order differential equation 

1 //" H"g' 2 g"H' n H'g'f\ M> ( f H'g'\ „ 
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in which primes denote derivatives with respect to g and z for the functions H(g), 
g{z) and f(z) respectively. With the confidence gained by the similarity between 
Eqs. (15) and (4), one can safely use the present treatment splitting Eq. (15) in 
two pieces 



W 2 (z) - 



W(z) 



V (z)-e, W = 



r 



Mf 



and 



AW 2 {z)- 



AW(z) 



M 



+ 2W(z) AW(z) = AV(z) - AE, AW 



H'g> 



MH 



(16) 



(17) 



which is similar to Eqs. (6)and (7), where À = e + AE and V e ff = Vq + AV. 

After all, it can be clearly seen that Eq. (16) is the one required for obtaining 
an explicit expression for W term used in Eq. (6) corresponding to an exactly 
solvable system considered in one-dimension (a — 7 = 0). However, to proceed 
further, the functions / and g should be solved as H, Q and are R known in 
principle. Now, equating like terms between the resulting expression in (15) and 
(14) gives 



where, from the definition of Q, 



f" M' f 

—L· _|_ M(E - V) 

f M f 1 



/(*) - (^) V2 



9' 



cxp 



1 f9(*) 



Q{g)dg 



(18) 
(19) 



Consideration of Eqs. (15) through (18) suggests a novel prescription 

AV{z)-AE = -ÇR[g{z)\ , 



(20) 



which, for plausible M and R functions, provides a reliable expression for g(z). 
It is remarked that in the constant mass case M — > 1 this procedure reduces 
to the well known formalism which has been throughly investigated [11] that, 
together with [9], justify our new proposal in solving PDMSE. The more detailed 
investigation of this treatment will be discussed elsewhere. 



5 Concluding Remarks 

In this Letter, a general method has been presented to address the question of 
corrections to the solution in one-dimension for a large class of iV— dimensional 
and exactly solvable PDMSE. We have also described how to extend the method 
to the case where the necessary function W(r) in (6) generating algebraically 
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solvable potentials in one dimension are present, which initiates calculations in 
the model leading to explicit expressions for the modificat ions due to both the 
use of physically plausible Zhu-Kroemer effective Hamiltonian (a = 7 = —1/2) 
and higher dimensional treatments. The main results are consistent with the 
other related works in the literature, which allow a non-perturbative treatment 
of these issues. 

Although, for clarity we have illustrated an application of the method for an 
easily accessible case of interest, it can be readily employed in various typical 
situations. In view of the importance in calculating such corrections in physics, 
we believe that the present model would serve as a useful toolbox to treat even 
more realistic situations which now occur in experimental observations with the 
advent of the quantum technology. 

References 

[1] Òzer O and Gònül B 2003 Mod.Phys. Lett. A18 2581; Gònül B 2004 Chmese 
Phys.lett. 21 1685; Gònül B and Koçak M 2005 Mod.Phys. Lett. A20 355; 
Gònül B, Çelik N and Olgar E 2005 Mod.Phys. Lett. A20 1683; Gònül B and 
Koçak M 2005 Mod.Phys. Lett. A20 1983; Gònül B, Kòksal K and Bakir E; 
quant-ph/0507098. 

[2] Gònül B 2004 Chmese Phys.lett. 21 2330. 

[3] Gònül B and Koçak M 2005 Chmese Phys.lett. 22 2742. 

[4] Quesne C 2005 Ann. Phys.(N. Y.) (to appear); quant-ph/0508216; Yu J, 
Dong S and Sun G 2004 Phys. Lett. A322 290; Chen G 2004 Phys. Lett. 
A329 22; Chen G and Chen D 2004 Phys. Lett. A331 312; Chen Z and 
Chen G 2005 Physica Scripta 72 11. 

[5] Gònül B, Gònül B, Tutcu D and Òzer O 2002 Mod.Phys. Lett. A17 2057; 
Gònül B, Òzer O, Gònül B and Üzgün F 2002 Mod.Phys. Lett. A17 2453. 

[6] BenDaniel D J and Duke C B 1966 Phys. Rev. 152 683; Bastard G 1981 
Phys. Rev. B24 5693; Zhu Q G, Kroemer H 1983 Phys. Rev. B27 3519; Li 
T, Kuhn K J 1993 Phys. Rev. 47 12760. 

[7] Gònül B, Òzer O, Koçak M, Tutcu D and Cançelik Y 2001 J. Phys. A34 
8271. 

[8] Plastino A R, Rigo A,Casas M, Garcias F and Plastino A 1999 Phys. Rev. 
C60 4318. 

[9] Bagchi B, Gorain P, Quesne C and Roychoudhury R 2005 Europhys. Lett. 
72 155. 



7 



[10] Abramowitz M and Stegun I A 1970 Handbook of Mathematical Functions 
(New York: Dover) 

[11] Bhattacharjie A and Sudarshan E C G 1962 Nuovo Cimento 25 864; Natan- 
zon G A 1979 Theor. Math. Phys. 38 146; Levai G 1989 J. Phys. A22 689; 
Levai G 1991 J. Phys. A24 131; Znojil M and Levai G 2001 J. Phys. A42 
1996; Bagchi B and Ganguly A 2003 J. Phys. A36 L161. 



8 



